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Transitions between inverse anticipatory, inverse complete and inverse lag synchronizations are 
shown to occur as a function of the coupling delay in unidirectionally coupled time-delay systems 
with inhibitory coupling. We have also shown that the same general asymptotic stability condition 
obtained using the Krasovskii-Lyapunov functional theory can be valid for the cases where (i) both 
the coefficients of the A{t) (error variable) and At — A{t — t) (error variable with delay) terms in 
the error equation corresponding to the synchronization manifold are time independent and (ii) the 
coefficient of the A term is time independent while that of the At term is time dependent. The 
existence of different kinds of synchronization are corroborated using similarity function, probability 
of synchronization and also from changes in the spectrum of Lyapunov exponents of the coupled 
time-delay systems. 
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Synchronization of chaotic oscillations is one 
of the most fundamental phenomena exhibited 
by coupled chaotic oscillators. Since the iden- 
tification of chaotic synchronization in identical 
systems, different kinds of synchronizations such 
as generalized, phase, lag, anticipatory and in- 
termittent synchronizations, and their variants, 
have been reported in the literature. Recent 
studies on synchronization have also been focused 
on the existence of inverse synchronization. In 
this connection, we have identified different kinds 
of inverse synchronization and transitions among 
them as a function of a single parameter in cou- 
pled time-delay systems with inhibitory coupling. 
It is now recognized that inhibitory couplings 
are important and widespread in nature, specif- 
ically in biological systems, neurophysiology and 
in many natural networks. We have also obtained 
a suitable stability condition for asymptotic sta- 
bility of the synchronized state using Krasovskii- 
Lyapunov functional theory. The results are cor- 
roborated using similarity function, probability of 
synchronization and from the changes in the spec- 
trum of Lyapunov exponents of specific coupled 
time-delay systems. 
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I. INTRODUCTION 



Synchronization is an interesting dynamical phe- 
nomenon exhibited by interacting oscillators in diverse 
areas of science and technology |l|, 0]. It has become 
an area of active research since the identification of syn- 
chronization in chaotic oscillators In recent years, 
different types of synchronization and generalizations 
have been reported both experimentally and theoreti- 
callyjllj 

Among them, inverse or anti-synchronization 
m, [Mm, [M [13, m, m constitutes an important class 
of synchronization, which is a typical feature of dynami- 
cal systems interacting through inhibitory coupling. 

By definition, inverse synchronization is the phe- 
nomenon 20, 21, 22, 23, 24, 25, 26, 21,28,,^] where the 
state vectors, say the drive x{t) and the response y{t), of 
the synchronized systems have the same absolute values 
but of opposite signs, that is x{t) — —y{t). While this is 
inverse complete synchronization, one may also identify 
inverse anticipatory synchronization where the response 
anticipates the drive, that is y{t) — —x{t + r), where 
T > is a constant. Similarly, when the response system 
lags the drive, one has inverse lag synchronization, that 
y{t) — ~x{t — r), T > 0. Just as in the case of regular 
synchronization, especially in nonlinear dynamical sys- 
tems with inhibitory coupling, one may encounter all the 
above types of inverse synchronization for appropriate 
ranges of the parameters. 
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The importance of inhibitory or repulsive coupling is 
well acknowledged in biological systems. It is a well es- 
tablished fact that couplin gs b etween neurons are both 
excitatory and inhibitory [30|. Ecological webs typi- 
cally have both positive and negative connections be- 
tween their components [3l|. Coupled lasers with neg- 
ative couplings have also been widely studied [32l|. The 
well known Swift-Hohenberg and Kuramoto-Sivashinsky 
equations include such a term [33[ . Currently, it has also 
been realized that a large class of natural networks also 
have inhibitory interactions among the nodes [s^, • 

The first experimental observation of inverse synchro- 
nization was demonstrated in coupled semiconductor 
laser diodes [20| , in which it was established that inverse 
synchronization was caused by nonresonant coupling be- 
tween the drive and the response lasers. It was also shown 
that switching between synchronization and inverse syn- 
chronization was possible by slightly changing the pump 
current of the drive laser [21| . Experimental observations 
and numerical simulations of synchronization and inverse 
synchronization of low frequency power drop-outs and 
jump-ups of chaotic semiconductor lasers were carried 
out in [231 . Inverse synchronization was also observed 
both experimentally and numerically in unidirectionally 
coupled laser systems with optical feedback [23| , as 
well as in a class of chaotic delayed neural networks [2^ 
and in coupled Ikeda systems with multi-feedback and 
multiple time-delays [25|. Inverse anticipating synchro- 
nization was demonstrated in coupled Ikeda systems [2^ . 
Further inverse retarted/lag synchronization and the role 
of parameter mismatch were discussed in (26j . However 
it may be noted that in none of the above studies the role 
of inhibitory coupling was investigated. 

Despite the fact that a vast amount of literature is 
now available on the phenomenon of synchronization, in- 
verse synchronizations have not been studied adequately, 
in particular with inhibitory couplings between interact- 
ing dynamical systems. In this paper, we report inverse 
synchronizations (inverse anticipatory, inverse complete, 
and inverse lag synchronizations) in unidirectionally cou- 
pled time-delay systems with inhibitory coupling. We 
also present a sufficient stability condition for asymp- 
totic stability of the synchronized state following the 
Krasovskii-Lyapunov functional approach for the cases 
where (i) the coefficients of the A{t) (error variable) and 
A,- = A(t — r) (error variable with delay) terms of the er- 
ror equation corresponding to the synchronization man- 
ifold are constant and (ii) the coefficient of the A term 
is constant and that of the A,- term is time dependent. 
We show that there is a transition from inverse antici- 
patory to inverse lag synchronization through complete 
inverse synchronization as a function of the delay time in 
the coupling. The tools to get these results are similar- 
ity function, probability of synchronization and largest 
Lyapunov exponents of the coupled time-delay systems. 

The plan of the paper is as follows. In Sec. II, we de- 
duce a sufficient condition for the asymptotic stability 
of the synchronized state for a system of unidirectionally 



coupled scalar delay differential equations with inhibitory 
coupling. We consider a piecewise linear delay differential 
equation as an example for the case where the coefficients 
of both the A and A.^ terms in the error equation are con- 
stant and demonstrate the existence of different types of 
inverse synchronization as a function of the coupling de- 
lay in Sec. III. Using the paradigmatic Ikeda system in 
Sec. IV for the case where the coefficient of the A,- term 
is time dependent while that of the other is time indepen- 
dent, we show that the same general stability condition is 
valid for the asymptotic stability of different types of in- 
verse synchronizations again as a function of the coupling 
delay. We also demonstrate these dynamical transitions 
through numerical analysis. Finally, in Sec. V, we sum- 
marize our results. 



II. COUPLED SYSTEM AND THE STABILITY 
CONDITION 

Consider the following unidirectionally coupled drive, 
x(t), and response, y{t), systems with inhibitory coupling 
of the form 

x{t) = ~ax{t) + bif{x{t-Ti)), (la) 
m - -ay{t) + b2f{y{t-n))-b3f{x{t^T2))ilh) 

where 61,62 and 63 are positive parameters, a > 0, ri 
and T2 are the feedback and the coupling delays, respec- 
tively. The nonlinear function f{x) is chosen to be a 
piecewise linear function which has been studied in de- 
tail recently [13, [H, [s?" 



X < -4/3 
-4/3 < a; < -0.8 
-0.8 < a; < 0.8 
0.8 < a; < 4/3 

X > 4/3, 




as the first example, and 



(2) 



/(^) 



sm a; 



(a;,-) = sin(a;(< — r)). 



(3) 



which is the well known Ikeda system [4(| , as the second 
example. 

Now the stability condition for the synchronization 
of the coupled time-delay systems, Eqs. ([1]), with the 
inhibitory delay coupling, —b3f{x{t — T2)), can be ob- 
tained as follows. The time evolution of the difference 
system (error function), associated with inverse synchro- 
nization, with the state variable A(t) = Xt2-ti + y(Oi 
where Xr2-Tx — x{t — {t2 — ti)), can be written for small 
values of A, by using the evolution Eqs. ((T|), as 



A = 



(4) 



= -aA + (61 - 62 - b3)f{x{t - T2)) + b2f{x{t - T2))A(^) 

The above evolution equation (O corresponding to the 
error function of the inverse synchronization manifold is 
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inhomogeneous and so it is difficult to analyse the system 
analytically. Nevertheless, the evolution equation can be 
written as a homogeneous equation 

A^-aA + b2f'ix{t-T2))Ar„ (6) 

for the specific choice of the parameters 

61 = 62 + 63- (7) 

Therefore, we will concentrate on this parametric choice. 
The inverse synchronization manifold A = Xt2-ti +2/ = 
corresponds to the following distinct cases: 

1. Inverse anticipatory synchronization occurs when 
T2 < Ti withy(t) = —x{t—T);T = T2—T1 < 0, where 
the state of the response system anticipates the in- 
verse state of the drive system in a synchronized 
manner with the anticipating time f (whereas in 
the case of direct anticipatory synchronization, the 
state of the response system anticipates exactly the 
state of the drive system, that is, y{t) = x{t — f)). 

2. Inverse complete synchronization results when T2 = 
Ti with y[t) = —x(t)]f = T2 — Ti =0, where the 
state of the response system evolves in a synchro- 
nized manner with the inverse state of the drive 
system (whereas in the case of complete synchro- 
nization, the state of the response system evolves 
exactly identical to the state of the drive system, 
that is, y{t) — x{t)). 

3. Inverse lag synchronization occurs when T2 > ti 
with y{t) = —x{t — f);f = T2 — Ti > 0, where the 
state of the response system lags the inverse state 
of the drive system in a synchronized manner with 
the lag time f (whereas in the case of direct lag 
synchronization, the state of the response system 
lags exactly the state of the drive system, that is, 
yit)=xit-T)). 

The synchronization manifold corresponding to Eq. ([6]) 
is locally attracting if the origin of the above error equa- 
tion is stable. Following the Krasovskii-Lyapunov func- 
tional approach [36. , .38] , we define a positive definite Lya- 
punov functional of the form 

V{t) = -A^+fi / A^{t + 9)de, (8) 

where fi is an arbitrary positive parameter, > 0. 

It is to be noted that for the above error equation 
the coefficient of the A term is always constant while that 
of the At- term can be time dependent. Hence one can 
obtain two cases depending on the choice of the nonlin- 
ear functional form f{x). If the derivative of the function 
turns out to be a constant as in the case of the piecewise 
linear function, Eq. ([2]), one obtains a constant coefficient 
for the At term. On the other hand, if the derivative 
still depends on time as in the case of the Ikeda system, 
Eq. ([3]), then the At- term always has a time dependent 



coefficient. In the following, we will show that the same 
general stability condition derived from the Krasovskii- 
Lyapunov approach can be valid for both the cases for the 
asymptotic stability of the different types of inverse syn- 
chronizations. However, there also arises an even more 
general situation where the coefficients of both the A and 
At- terms are time dependent. In this case the arbitrary 
positive parameter fj, in the Lyapunov functional is no 
longer a positive constant. We have also designed a suit- 
able coupling form for this rather general situation, for 
which we can show that the same general stability con- 
dition arrived from the Krasovskii-Lyapunov approach is 
still valid for the asymptotic stability of the synchronized 
state; these results will be published elsewhere. 

Note that from Eq. ([8]), V{t) approaches zero as A — * 
0. Hence, the required solution A = to the error equa- 
tion, Eq. is stable only when the derivative of the 
Lyapunov functional V{t) along the trajectory of Eq. ^ 
is negative. This requirement results in the condition for 
stability as 

r(/i) = A^l{a - ^) > bjf'ixiit - T2)f. (9) 

Again r(/i) as a function of fi for a given f'{x) has an 
absolute minimum at ^ — a/2 with Tmin — ■ Since 
r > r„ii„ = a^, from the inequality ([5]), it turns out that 
a sufficient condition for asymptotic stability is 

a>\h2f(x{t-T2))\. (10) 

This general stability condition indeed corresponds to the 
stability condition for inverse anticipatory, complete in- 
verse and inverse lag synchronizations for suitable values 
of the coupling delay T2 corresponding to a fixed value 
of the feedback delay t\ for both the piecewise linear 
and Ikeda time-delay systems corresponding to the cases 
where the coefficient of the A,- term in the error equation 
is time independent and time dependent, respectively. 

Further, it is interesting to note that if one substitutes 
?/ — > y = —y in Eq. ([ija), then the coupling becomes exci- 
tatory for the choice of functional forms we have chosen. 
This is exactly the case we have studied in 38] , where di- 
rect anticipatory, complete and lag synchronizations exist 
as a function of the coupling delay. However, one cannot 
obtain inverse (anticipatory, complete and lag) synchro- 
nization with excitatory coupling or direct (anticipatory, 
complete and lag) synchronization with inhibitory cou- 
pling for the chosen form of the unidirectional nonlinear 
coupling because of the nature of the parametric relation 
between h\ , 62 and 63 and the stability condition pU|) . 

Now from the form of the function /(x) in Eq. ^ for 
the piecewise linear time-delay system, one can obtain a 
less stringent stability condition as [3^ 

a>62, (11) 

while 

a>1.5&2, (12) 
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FIG. 1: (a) The first ten maximal Lyapunov exponents Xmax 
of tlie piecewise linear time-delay system ((H^) and ([2}) for 
the parameter values a — 1.0, b = 1.2, t £ (2,29) and (b) 
Hyperchaotic attractor for the delay time t = 8.0 with two 
positive Lyapunov exponents for the above values of the other 
parameters. Here n — r. 
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FIG. 2; (Color online) The time trajectory of the variables 
x{t),y{t) and -yit) for a = 1.0, foi = 1.2, 62 = 0.6,63 = 
0.6, Ti = 8.0 and T2 — 6.0 of the coupled piecewise linear 
time-delay systems indicating IAS. 



is the most general condition specified by ([TO]) for asymp- 
totic stability of the synchronized state A = 0. 

Correspondingly, one can obtain the stability condition 
for the coupled Ikeda systems as 



a > I62 cos{x{t — T2))| 



(13) 



III. INVERSE SYNCHRONIZATIONS IN THE 
COUPLED PIECEWISE LINEAR TIME-DELAY 
SYSTEMS ((P AND ^) 

Here, we will show the existence of inverse anticipatory, 
inverse complete and inverse lag synchronizations as a 
function of the coupling delay in the coupled piecewise 
linear time-delay systems (IT])-(I2|) with inhibitory coupling 
as an illustration for the case where the coefficients of 
both the A and A^ terms in the error equation ^ are 
constant. 



A. Inverse anticipatory synchronization, IAS 

In this section, we present the existence of inverse an- 
ticipatory synchronization (IAS) for the values of the cou- 



pling delay T2 < ri, and for fixed values of the other 
system parameters. In particular, we have fixed the pa- 
rameters as a = 1.0, 5i = 1.2, ri = 8.0, T2 ~ 6.0 and the 
other two parameters 62 and 63 arc fixed according to 
the parametric condition ([7]) . The first ten maximal Lya- 
punov exponents Xmax of the uncoupled piecewise linear 
time-delay system in the range of delay time r G (2, 29) 
for the above choice of parameters are shown in Fig. [IJi. 
It is clear from Fig. [1^ that for t > 0.5 at least two of 
the Lyapunov exponents are positive and that the sys- 
tem is hyperchoatic. As an illustration, the hyperchaotic 
attractor with two positive Lyapunov exponents for the 
value of the delay time n = r = 8.0 is plotted in Fig.[T]D. 
The time trajectories of the variables x{t), y{t) and —y{t) 
of the coupled piecewise linear time-delay systems (II|)-(I1]) 
are plotted in Fig. [2] depicting the existence of IAS for 
62 = 0.6, for which the more general stability condition, 
a > 1.5&2, is satisfied. It is evident from the figure that 
the state of the response system y{t) anticipates the in- 
verse state of the drive system x{t). This is made visually 
clear by plotting the inverse state of the response, that 
is —y{t). From Fig. [21 it is clear that the inverse state 
of the response {—y{t)) anticipates the state of the drive 
x{t), thereby illustrating the existence of IAS between 
the drive and the response systems. 

The existence of IAS for the above choice of param- 
eters is further confirmed using the notion of similarity 
function S{b2), probability of synchronization $(^2) and 
also from the spectrum of Lyapunov exponents of the 
coupled time-delay systems. Now we use the notion of 
similarity function, introduced in [llj for characterizing 
the lag synchronization, to characterize the existence of 
IAS. The similarity function, Sa(b2), is defined as a time 
averaged difference between the variables x{t) and —y{t) 
(with mean values being subtracted) taken with the time 
shift T, 



{[y{t-r)+xm 
[{x^t)){y^tW/^ 



where (x) means time average over the variable x. If the 
signals x{t) and —y{t) are independent, the difference be- 
tween them is of the same order as the signals themselves. 
If x{t) = y{t), as in the case of complete synchronization, 
the similarity function reaches a minimum, Sa(b2) = 0, 
for r = 0. On the other hand, if 6*0 (^2) = for the case 
T 7^ 0, there exists a time shift r between the two signals 
x{t) and —y(t) such that y{t — r) = —x(t), demonstrat- 
ing inverse anticipatory synchronization. The similarity 
function, 6*0 (&2), as a function of the parameter &2 is 
shown in Fig. [3^, which clearly indicates that the value of 
Sa{b2) oscillates with finite amplitude (5a (^2) > 0) above 
the value of the control parameter &2 ~ 0.97, indicating 
the desynchronized evolution between the response, y{t), 
and the drive, x{t), systems. For the value of the control 
parameter 62 < 0.97, the similarity function become zero 
iSaib2) — 0) indicating that the response system antici- 
pates the inverse state of the drive system confirming the 
existence of IAS. 
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FIG. 3: (a) Similarity function, Sa{b2), (b) Probability of 
synchronization, ^{b2) and (c) First eight maximal Lyapunov 
exponents, Xmax, of the coupled piecewise linear time-delay 
systems as a function of the control parameter 62, indicating 
the existence of IAS for 62 < 0.97. 



It is to be noted that one can observe exact inverse 
anticipatory synchronization in the present study even 
for the choice of the parameter &2 for which only the less 
stringent stability condition pT|) is satisfied in contrast 
to our earlier studies [38| , where one can observe only ap- 
proximate direct anticipatory synchronization under this 
condition while exact synchronizations are observed for 
the values of the parameters satisfying the general stabil- 
ity condition for the coupled piecewise linear time-delay 
systems. This is due to the fact that in the present study, 
we have chosen the value of the delay time in the drive 
system, x{t), as r — 8, for which the system has only 
two positive Lyapunov exponents (as may be observed 
from the Fig. [T]). This implies that the number of trans- 
versely unstable manifolds is less and hence they are all 
stabilized even for the least values of the parameters sat- 
isfying even the less stringent stability condition and so 
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FIG. 4: (Color online) The time trajectory of the variables 
x(t),y(t) and ~y{t) of the coupled piecewise linear time-delay 
systems indicating ICS for the coupling delay r2 = n = 8.0, 
while the other parameter values are the same as in Fig. (2] 



one can observe exact inverse synchronization. This is in 
contrast to our earlier studies on direct anticipatory syn- 
chronization [sll, where we have chosen r = 25, for which 
the drive system exhibits more than seven positive Lya- 
punov exponents. Correspondingly there exists a large 
number of transversely unstable manifolds and hence the 
general stability condition is required to be satisfied in 
order to stabilize all the transversely unstable manifolds 
to obtain exact synchronization. 

As inverse anticipatory synchronization is a special 
case of generalized synchronization, the existence of IAS 
can also be further confirmed using the auxiliary sys- 
tem approach by augmenting the coupled piecewise lin- 
ear time-delay systems (([1]) and ©) with an additional 
auxiliary system for the variable z{t) identical to the re- 
sponse system, satisfying the equation 

Z{t) = -az{t) + h2f{z{t - n)) - b^f{x{t - T2)). (15) 

Now, the existence of IAS can be characterized by us- 
ing the probability of synchronization 38], $(62), calcu- 
lated between the response, y{t), and the auxiliary sys- 
tems, z(i), which can be defined as the fraction of time 
during which \y{t) — z{t)\ < e occurs, where e is a small 
but arbitrary threshold. The probability of synchroniza- 
tion, $(62), remains zero for the value of the control pa- 
rameter 62 > 0.97 as shown in Fig. [3]d, where there is no 
synchronization between the response, y{t), and the aux- 
iliary system, z{t). On the other hand, for 62 < 0.97, the 
probability of synchronization attains the value of unity 
for the chosen threshold value for e, clearly indicating 
the existence of complete synchronization between the 
response, y{t), and the auxiliary system, z{t) (we have 
fixed the threshold value at e = 10"^" throughout the 
manuscript). Correspondingly, there exits IAS between 
the coupled drive, x{t), and the response, y(i), systems 
with inhibitory coupling. 

Further, the existence of IAS can also be characterized 
by the changes in the spectrum of the Lyapunov expo- 
nents of the coupled systems, ([I])-©. The first eight 
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FIG. 5: (a) Similarity function, S'c(&2), (b) Probability of 
synchronization, "1>(62) and (c) First eight maximal Lyapunov 
exponents, Xmax, of the coupled piecewise linear time-delay 
systems as a function of the control parameter 62, indicating 
the existence of ICS for 62 < 0.97 and for the coupling delay 
T2 = n = 8.0. 



largest Lyapunov exponents of the coupled piecewise lin- 
ear time-delay system is shown in Fig.[3J:. The two largest 
Lyapunov exponents of the drive system, x{t), remain 
positive, while the ones corresponding to the response 
system, y(t), decrease in their value as a function of the 
parameter 62. The least positive Lyapunov exponent of 
the response system becomes negative at 62 = 1.06, while 
the largest positive Lyapunov exponent of the response 
system becomes negative at the value of 62 = 0.97 con- 
firming the onset of IAS at 62 — 0.97. 



B. Inverse Complete synchronization, ICS 

The synchronization manifold, A = Xt-2-ti + y, be- 
comes an ICS manifold for ti = T2 — 8.0. The time tra- 




2500 



2520 



2540 



2560 



2580 



2600 



FIG. 6: (Color online) The time trajectory of the variables 
x(t),y(t) and ~y{t) of the coupled piecewise linear time-delay 
systems indicating ILS for the value of the coupling delay 
T2 = 10.0, while the other parameter values are the same as 
in Fig. m 



jectory of the variables x{t),y{t) and —y{t) are plotted 
in Fig. 21 depicting the existence of ICS for the same val- 
ues of the other parameters as in Fig. [2l It is clear from 
the Fig. [5] that the state of the response system, y{t), 
evolves in synchrony with the inverse state of the drive 
system x{t). For the ease of clear visualization of this 
phenomenon, the inverse of the response —y{t) is shown 
as filled circles in Fig. [2] and hence the inverse state of 
the response —y(t) clearly evolves in synchrony with the 
state of the drive x{t), there by illustrating the existence 
of ICS between the drive, y{t), and the response, x{t), 
systems. 

We have also calculated the similarity function for the 
ICS defined as 



{[y{t) + x{m 

[{xHt)){yHtW/^-' 



(16) 



as a function of the parameter 62 in Fig. The similar- 
ity function, Sc{b2), oscillates with a finite value above 
zero for the value of the control parameter 62 > 0.97 as 
shown in Fig. [5|i, where there does not exist any cor- 
relation between the interacting systems. However, it is 
evident from Fig.Oa, that the similarity function becomes 
zero for the values of 62 < 0.97 indicating the existence 
of exact ICS in the coupled piecewise linear time-delay 
system even for the value of the parameter satisfying the 
less stringent stability condition as discussed in the pre- 
vious section. We have also characterized the existence 
of ICS using the probability of synchronization ($(62)) 
as shown in Fig. [Sb, which clearly shows that the prob- 
ability of synchronization becomes unity for 62 < 0.97 
depicting the existence of ICS between the coupled drive, 
x{t), and the response, y{t), systems with inhibitory cou- 
pling. However, for 62 > 0.97, the value of the probabil- 
ity of synchronization becomes zero, $(62) = 0, indicat- 
ing that there exist no correlation between the coupled 
systems for this range of parameters. The eight largest 
Lyapunov exponents of the coupled time-delay systems is 
shown in Fig.[5t. It shows that the two largest Lyapunov 
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FIG. 7: (a) Similarity function, £'((62), (b) Probability of syn- 
chronization, $(62) and (c) First eight maximal Lyapunov 
exponents, Xmax, of the coupled piecewise linear time-delay 
systems as a function of the control parameter 62, indicating 
the existence of ILS for 62 < 0.97 and for the coupling delay 
T2 = 10.0. 



exponents of the drive system remain positive while the 
largest positive Lyapunov exponent of the response sys- 
tem attains a negative value at 62 = 0.97, confirming the 
existence of ICS for 62 < 0.97. 



FIG. 8: (a) The first eleven maximal Lyapunov exponents 
Xmax of the Ikeda time-delay system ((H^) and for a = 
1.0,6 = 5 r G (2,25) and (b) Hyperchaotic attractor for the 
delay time ti = t — 4.0 with three positive Lyapunov expo- 
nents for the above values of the other parameters. 



system y{t) lags the inverse state of the drive system 
x{t). In order to made this clear visually, the inverse 
state of the response —y{t) is plotted in Fig. [5] and hence 
the inverse state of the response — lags the state of 
the drive x(t), there by illustrating the existence of ILS 
between the drive and the response systems. 
The similarity function for ILS defined as 

02., ^ {[V{t + T)+X{t)f) 

is plotted in Fig. [7^ as a function of &2 ■ The similarity 
function oscillates with a finite amplitude for &2 > 0.97 
as in the other cases, as there is no synchronous evolu- 
tion among the coupled systems. However, the similarity 
function, Si{b2), reaches zero for 62 < 0.97, confirming 
the existence of ILS in the coupled time-delay systems. 
Similarly, the value of the probability of synchroniza- 
tion, $(62), calculated between the response, y{t), and 
the auxiliary, z(t), systems attains unity for 62 < 0-97, 
confirming the existence of complete synchronization be- 
tween them as shown in Fig. [7|d. Correspondingly, there 
exists ILS between the drive, x(t), and the response, y{t), 
systems. The two largest Lyapunov exponents of the 
drive system remain positive as shown in Fig. [7b, while 
the largest positive Lyapunov exponent of the response 
system becomes negative at 62 = 0.97, confirming the 
existence of ILS between the drive and the response sys- 
tems. 



C. Inverse lag synchronization, ILS 

Now, we will demonstrate the existence of ILS in the 
coupled piecewise linear time-delay system, ([l])-©, for 
the value of coupling delay T2 > ti and for the same 
values of the other parameters as in Sec. IIII Al The time 
series of the variables x{t), y{t) and —y{t) of the coupled 
systems are shown in Fig. [5] for &2 = 0.6 and = 10.0, 
depicting the existence of ILS in the coupled systems. It 
is evident from the figure that the state of the response 



IV. INVERSE SYNCHRONIZATIONS IN THE 
COUPLED IKEDA TIME-DELAY SYSTEMS ({iD 
AND dSD) 

In this section, as an illustration for the case where 
the coefficient of the A term in the error equation ^ is 
constant, while that of the A,- term is time-dependent, 
we will demonstrate the existence of inverse anticipatory, 
inverse complete and inverse lag synchronizations as a 
function of the coupling delay in the coupled Ikeda time- 
delay systems ([1]) and Q with inhibitory coupling. 
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FIG. 9; (Color online) The time trajectory of the variables 
x{t),y{t) and ~y{t) for a = 1.0, 6i = 5.0,^2 = 2.0,63 = 
3.0, Ti = 4.0 and T2 = 3.0 of the coupled Ikeda time-delay 
systems indicating IAS. 

A. Inverse anticipatory synchronization, IAS 

We have fixed the values of the parameters of the 
coupled Ikeda time-delay systems (([1]) and ([3])) as a = 
1.0,61 = 5,ri — 4.0 and T2 = 3.0, while the other two 
parameters 62 and 63 are fixed according to the paramet- 
ric condition bi =62 + 63. To appreciate the chaotic and 
hyperchaotic nature of the uncoupled system, we present 
in Fig. [8^ the first eleven largest Lyapunov exponents for 
the above values of the parameters in the range of delay 
time T € (2, 25) where several of them take positive val- 
ues. As an illustration, the hyperchaotic attractor with 
three positive Lyapunov exponents of the uncoupled sys- 
tem for Ti = T = 4.0 is depicted in Fig. [Hd. 

Now, wc will demonstrate the existence of IAS for the 
value of the coupling delay T2 less than that of the feed- 
back delay ti and for the values of the other parameters 
satisfying the stability condition (|13p . It is evident from 
Fig- that the maximum value of x(t) does not exceed 
Xmax = 5. Correspondingly x(t - T2)max = 5. As a 
consequence the stability condition (|13p can be written 
as 

a > 62 cos(5) = 0.2 8 462. (18) 

Then, one can obtain an asymptotically stable synchro- 
nized state for the values of parameters satisfying the 
above stability condition. The time trajectories of the 
variables x{t),y(t) and —y{t) depicting the existence of 
IAS are plotted in Fig. [S] for the value of the parame- 
ter 62 = 2.0 satisfying the stability condition (fT5|) . The 
minimum of the similarity function, 6*0(62), defined by 
Eq. becomes zero for 62 < 2.88 as shown in Fig. [TUb. 
indicating the existence of IAS in the coupled Ikeda time- 
delay system. 

The existence of IAS is further characterized by the 
probability of synchronization, ^(b2), for complete syn- 
chronization between the response, y(t), and the auxil- 
iary, z{t), systems by augmenting the coupled Ikeda time- 
delay systems (([T} and ([3])) with an additional auxiliary 
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FIG. 10: (a) Similarity function, Si{b2), (b) Probability of 
synchronization, "1>(&2) and (c) First nine maximal Lyapunov 
exponents, Xmax, of the coupled Ikeda time-delay systems as a 
function of the control parameter 62, indicating the existence 
of IAS for 62 < 2.88 and for the coupling delay T2 = 3.0. 
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FIG. 11: (Color online) The time trajectory of the variables 
x{t),y{t) and -yit) for a = 1.0, foi = 5.0, 62 = 2.0,63 = 
3.0, Ti = 4.0 and t2 = 4.0 of the coupled Ikeda time-delay 
systems indicating ICS. 
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FIG. 12: (a) Similarity function, Si{b2), (b) Probability of 
synchronization, "I>(&2) and (c) First nine maximal Lyapunov 
exponents, Xmax , of the coupled Ikeda time-delay systems as a 
function of the control parameter &2, indicating the existence 
of ICS for b2 < 2.88 and for the coupling delay T2 = n = 4.0. 




FIG. 13: (Color online) The time trajectory of the variables 
x{t),y{t) and -y{t) for a = 1.0, 6i = 5.0,^2 = 2.0,63 = 
3.0, Ti = 4.0 and T2 = 5.0 of the coupled Ikeda time-delay 
systems indicating ILS. 



B. Inverse complete synchronization, ICS 

The time series of the variables x{t),y{t) and —y{t) 
indicating the existence of ICS, for the value of the cou- 
pling delay T2 = ti = 4.0, and for the same values of the 
other parameters satisfying the stability condition (fT5)) 
as in the previous section, are shown in Fig. 1111 The 
similarity function, Sc{b2), given by (jl6p . for ICS shown 
in Fig. [T^ indicates that the minimum of Sc{h2) = 
for 62 < 2.88 depicting the existence of ICS in the cor- 
responding range of 62. It is also evident from the prob- 
ability of synchronization (Fig. WBp ] that $(62) = 1 for 
62 < 2.88 confirming the existence of ICS between the 
drive and the response systems. This transition from a 
desynchronized state to ICS for the values of parameters 
satisfying the stability condition (|18p is also confirmed 
from the changes in the spectrum of the Lyapunov expo- 
nents of the coupled Ikeda systems as shown in Fig. [I2b . 
The largest positive Lyapunov exponents of the response 
system become negative at 62 = 2.88 confirming the exis- 
tence of ICS for 62 < 2.88, while the largest three positive 
Lyapunov exponents of the drive system remain positive. 



system for the variable z{t) identical to the response sys- 
tem as in Eq. (|15p . Figure. [TUb indicates that the value of 
the probability of synchronization becomes <&(62) = 1-0 
for &2 < 2.88, confirming the existence of IAS between 
the coupled drive, xli:), and the response, ?/(i), systems. 
The existence of IAS for 62 < 2.88 is further confirmed 
from the changes in the spectrum of the largest Lyapunov 
exponents of the coupled Ikeda time-delay systems (([T]) 
and (l3|)). The nine largest Lyapunov exponents, Xmax, is 
shown in Fig. [TUb as a function of 62- While the largest 
three positive Lyapunov exponents of the drive system 
remain positive, the largest positive exponents of the re- 
sponse system gradually become negative as a function 
of 62 and finally they become negative at 62 — 2.88, con- 
firming the existence of IAS for 62 < 2.88 satisfying the 
stability condition 



C. Inverse lag synchronization, ILS 

For the coupling delay T2 = 5.0 greater than the feed- 
back delay ti — 4.0, the synchronization manifold A = 
corresponds to the ILS manifold. The time trajectories 
of the variables x{t), y{t) and —y{t) of the coupled Ikeda 
time-delay systems (Fig. [T5)) clearly depict the ILS for 
the above choice of the coupling delay (the other pa- 
rameters are fixed as in Sec. IIV Al satisfying the stability 
condition (jlSp ). The minimum of the similarity func- 
tion for ILS turns out to be 5; (62) = as shown in 
Fig. [Mb for 62 < 2.88 indicating the existence of ILS 
for 62 < 2.88. Similarly, the value of the probability 
of synchronization (Fig. 114b ) becomes unity in the cor- 
responding range of 62 confirming the existence of ILS 
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FIG. 14: (a) Similarity function, Si{b2), (b) Probability of 
synchronization, "I>(&2) and (c) First nine maximal Lyapunov 
exponents, Xmax , of the coupled Ikeda time-delay systems as a 
function of the control parameter &2, indicating the existence 
of ILS for &2 < 2.88 and for the coupling delay T2 = 5.0. 



between the coupled Ikeda time-delay systems. The ex- 
istence of ILS for 62 < 2.88 is also further confirmed 
by the changes in the spectrum of the Lyapunov expo- 
nents of the coupled Ikeda time-delay systems as shown in 
Fig. [Hb . The largest three positive Lyapunov exponents 



of the drive system remain unchanged, while that of the 
response system decrease in their values as a function of 
62 and they become negative at 62 = 2.88 confirming the 
emergence of ILS in coupled Ikeda time-delay systems. 



V. SUMMARY AND CONCLUSION 

In this paper, we have shown the transition from in- 
verse anticipatory to inverse lag via inverse complete syn- 
chronization as a function of the coupling delay T2 for 
fixed value of the feedback delay ri in unidirectionally 
coupled time-delay systems with inhibitory coupling. We 
have also arrived a suitable stability condition for the 
asymptotic stability of the synchronized states using the 
Krasovskii-Lyapunov functional theory. We have demon- 
strated that the same general stability condition result- 
ing from Krasovskii-Lyapunov functional approach can 
be valid for two different cases, where (i) both the co- 
efficients of the A and A,- terms of the error equation 
corresponding to synchronization manifold are constants 
and (ii) the coefficient of the A,- term is time dependent 
while that of the other is time independent using suitable 
examples. The existence of different types of inverse syn- 
chronizations are corroborated using similarity function, 
probability of synchronization and from the changes in 
the spectrum of the largest Lyapunov exponents of the 
coupled time-delay systems. We have also designed suit- 
able couplings for the case where both coefficients of the 
A and A^ terms are time dependent to show the validity 
of the same general stability condition (jlOp resulting from 
the Krasovskii-Lyapunov functional theory, the results of 
that will be published in a forthcoming paper. 
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